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ABSTRACT 

We measure the intrinsic alignments of the blue and the red galaxies sepa- 
rately by analyzing the spectroscopic data of the Sloan Digital Sky Survey Data 
Release 6 (SDSS DR6). For both samples of the red and the blue galaxies with 
axial ratios of b/a < 0.8, we detect a 3a signal of the ellipticity correlation in 
the redshift range of < z < 0.4 for r-band absolute (model) magnitude cut of 
Mr < —19.2 (no K correction). We note a difference in the strength and the 
distance scale for the red and the blue galaxy correlation r]2D{'^)'- For the bright 
blue galaxies, it behaves as a quadratic scaling of the linear density correlation 
of ^(r) as rj2D{j') oc i'^{r) with strong signal detected only at small distance bin 
of r < 3/i^^Mpc. While for the bright red galaxies it follows a linear scaling 
as ri2D{r) oc ^(r) with signals detected at larger distance out to r ~ Gh^^Mpc. 
We also test whether the detected correlation signal is intrinsic or spurious by 
quantifying the systematic error and find that the effect of the systematic error 
on the ellipticity correlation is negligible. It is finally concluded that our results 
will be useful for the weak lensing measurements as well as the understanding of 
the large scale structure formation. 

Subject headings: cosmology:observations — methods:statistical — large-scale 
structure of universe 

It has been long suspected that the galaxies may have preferential directions in the 
orientations of their elliptical shapes (or spin axes) with being locally correlated between the 
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neighbors, which is often called the galaxy intrinsic alignments. Sir ice the galaxy intr insic 



alignments could in principle contaminate weak lensing signals (e.g.. iHirata et al.ll2007l . and 
references therein), the measurement of its existence and strength from observations has 
been a prior interest especially of the weak lensing community. 



It has been claimed by iMandelbaum et al.l (120061 ) that no signal of ellipticity correla- 
tion between the neighbor galaxies was det ected in the spectroscopic data from the Sloan 
Digital Sky Survey (SDSS. lYork et al.ll2000l ). although they mentioned that their results are 
statistically consisten t with previous ob s ervational evidences f or the existence of the galaxy 
intrinsic alignments (IBrown et al.ll2002l : iHeymans et al.ll2004l ). Extrapolating their results 
to redshift z = 1 relevant for the weak lensing shear, they have drawn a tentative conclusion 
that the intrinsic ellipticity correlations between the neighbor galaxies would not cause any 
significant errors for the estimation of the power spectrum amplitude of the initial density 
field. 

Very recently, however, a numerical analysis based on the high-resolution N-body simula- 
tion has demonstrated with high statistical power that strong ell ipticity correlat ions between 
the dark matter halos exist at large distances out to 10/i~^Mpc (ILee et al.ll2007l ). Therefore, 
it is necessary and important to reexamine observationally the existence and significance of 
the galaxy intrinsic alignments by analyzing the SDSS data. 

The physical mechanism for the generation of the galaxy intrinsic alignments differ be- 
tween the blue and the red galaxies. For the blue spiral galaxies which are usually rotationally 
supported, the intrinsic alignments refer to the correlations of the spin axes originated from 
the initial tidal field. The linear tidal torque theory explains that the intrinsic correlation 
should exist only locally between close galaxy pairs, whose fun ctional form is p redicted to 
be a quadratic scaling of the linear density correlation function (jPen et al.ll2000l ). 



While for the bright red galaxies which are not flattened by rotation, the intrinsic align- 
ments refer to the correlations between the major principal axes of the galaxy's elliptical 
shapes which are determined in observations in terms of the anisotropy in their stellar distri- 
bution. It has yet to be fully understood what the origin of the anisotropic steller distribution 
of the red galaxies is and how it is related to the shapes of the host halos. But the standard 
model based on the cosmic web theory suggests that it may be due to the anisotropic accre- 
tion and infall of materials and gas into the host halo along the halo's major pr incipal axes 



which are in turn elongated with t he large scale filamentary matter distribution (jWestlll994 



Bond. Kofman. fc Pogosyaru Il996l ). Given this anisotropic infall scenario, the red galaxies 
may have built-in memory of the large-scale filamentary structures and thus their ellipticities 
might be correlated on larger scales than the spin axes of the blue galaxies. 
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Even for the blue galaxies, it was claimed that the growth of th e non-Gaussianity in the 



density field would cause large-scale correlations of their spin axes (iHui fc Zhang||2002l ). To 
account for the large-scale intrinsic correlations between the neighbor galaxies, we use the 
following formula: 

where r is the three dimensional separation distance of a galaxy pair and ^(r; R) is the two- 
point correlation function of the linear density field smoothed on the Lagrangian galactic 
scale R. Here, the two correlation parameters a\ and ^ni represent the strength of the small- 
scale and the large-scale correlation, respectively. If ^ni = and a\ > 0, the galaxy intrinsic 
correlation will behave as a quadratic scaling of the density correlation function so that 
Vsoi'r) should rapidly diminishes to zero as r increases, as the linear tidal torque theory 
predicts. In contrast, if e„i > and ai = 0, the galaxy's intrinsic correlation will follow a 
linear scaling of ^(r) so that rjsnir) should be non-negligible even at large distance. 



In fact, we have suggested this formula (eq. [T]) in our previous work (ILee &: PenI 120071 ) 
to account for the non-Gaussianity effect on the galaxy's spin-spin correlations which are 
tidally induced. Here we use the same analytic model to describe the intrinsic ellipticity- 
ellipticity correlations of the galaxies, assuming that since the anisotropic infall of materials 
that resulted in the ellipticities of the galaxies are also nonlinear manifestation of the tidal 
field and thus can be approached with the same analytic techniques. Throughout this Letter, 
the galaxy intrinsic alignments refer to the intrinsic ellipticity-ellipticity alignments of the 
galaxies, unless otherwise stated. 

In practice, what can be measured is the intrinsic alignments not of the galaxy's three 
dimensional shapes but of its two dimensional shapes projected onto the plane of the sky. 



According to lLee &: Peru (1200 ll ). if the two dimensional projection effect is taken in to account 
properly, the correlation parameter has to be multiplied by a factor of 5/4. Thus, the two 
dimensional projected intrinsic correlation function is modified from equation ([1]) simply as 

25 , e{r;R) , 5 gr; R) 

Note here that although ?72_d('") represents the intrinsic correlations of the galaxy's two di- 
mensional projected ellipticities, the distance r represents the three dimensional separation. 
It is also worth mentioning here that for the blue galaxies the intrinsic correlations of the 
galaxy's two dimensional projected ellipticities are in fact same as that of the galaxy's two 
dimensional projected spin axes since the spin axes are believed to be orthogonal to the 
galaxy's major axes. 



To measure 'r]2D{r) from the real universe, we use a dataset of galaxies at redshift 
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< z < 0.4 downloaded from the SDSS DR6 website ( http:/ /www.sdss.org/dr67| . A 
total of 640647 galaxies are found to be in this redshift range. For each galaxy, we obtain 
information on right ascension (a), declination (5), redshift (z), position angle (p), major-to- 
minor axis ratios {b/a), the u — g and g — r colors, and r-band model magnitude {Mr). Among 
640647 galaxies, the information on p and b/a are available only for 639727 galaxies, from 
which we also select only those galaxies with b/a < 0.8 since for the nearly face-on galaxies 
with b/a > 0.8 the measurements of the direction of their major axes may suffer from large 
uncertainty. A total of 434849 galaxies are found to satisfy the criterion of b/a < 0.8. 



Using the empirical u — r color separator suggested by IStrateva et al.l (l200ll ) for the 
SDSS data, we classify the galaxies into the red (m — r > 2.22) and the blue {u — r < 2.22) 
sample. Then, we select only bright galaxies by applying the absolute magnitude cut of 
Mr < —19.2. Since the SDSS calibration uses the inverse hyperbolic sine magnitudes (asinh) 
(ILupton et al.lll999l ). we calculate Mr as the asinh magnitudes using the information on the 
softening parameter given in the SDSS web site. Finally, we end up having two SDSS samples 
which consist of a total of 87188 blue and 283972 red galaxies brighter than this magnitude 
cut of Mr < -19.2. 

For each galaxy, we determine the direction of the major axis, e, in the equatorial 
coordinate system from the given information on p, a and 6: Let us consider two galaxies in 
a given sample whose major axis directions are found as and e^. First, we calculate the unit 
projected separation vector as d = d/|d| with d = f j — fj, where f j and fj represent the unit 
position vectors of the i-th and the j-th galaxy, respectively. The separation vectors projected 
onto the plane of the sky at the position of the two galaxies are given as dj = d — (d • rj)rj 
and dj = d — (d ■ respectively. Then, the intrinsic correlations of the two dimensional 

projected ellipticities between the galaxies with separation of r is calculated as 



V2D[r) 



^cos^ (7i-7i) - J, 



with 



7j = tan 



fr,: X di) ■ e,; 



(3) 



Here, the angle represents the direction of projected onto the plane of sky at the position 
of Tj. Note that for a galaxy pair with small distance, the quantity A7 = 7j — 7^- is actually 
identical to the difference in the position angle Ap. To determine the three dimensional 
separation distance, r, of a galaxy pair, we assume a flat A-dominated cosmology with matter 
density Qm = 0.25, va cuum energy density Q\ = 0.75 and Hubble constant Hq = 100/ikm 
s-^Mpc^i with h = l Jeianton et al.ll2003h . 



Then, we measure ri2D{r) for the selected red and the blue galaxies, separately. Then, 
we compare the observational results with the analytic model ([2]), determining the best- 
fit values of the two parameters ai and ^ni with the help of the x^-minimization. In the 
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fitting pr ocedure, the values of the two parameters are confined to0<ai<0.6 and < 
^ni < O. KjLee fc Pen! l200ll . 120071 ) . The unmarginalized e r rors i n the fitting parameters are 
calculated by the formula given in lBevington fc RobinsonI (119961 ): The error in the parameter 
a\ is found as aa = Aai^y2{xi — 2Xmm + xD- Here Aa\ represents a fitting step-size, Xmm 



the minimum value of x^y Xi 



X (c^io + Aai) and Xi = X {o-io — Aai). where aio is defined 
Xmin = x'^{(^io)- The unmarginalized error cXg in the parameter ^ni is also calculated in a 
similar manner. 



F or the analytic model (eq.[2]), we use the ACDM power spectrum given by lBardeen et al, 
(Il986l ). For the Lagrangian scale R of each galaxy, we first calculate the luminosity L of 



each galaxy from its absolute r- band magnitude and th en convert it to the galaxy mass scale 
using the relation of L ~ M°'^^ (jVale &: Ostrikerll2006l ). Then, the Lagrangian galactic scale 
R is found as a top-hat radius that encloses the mean mass M averaged over the galaxies of 
each sample. When we calculate the absolute magnitude of each galaxy, no K-correction is 
considered: If our goal were to investigate how ri2D{f') changes in narrow absolute magnitude 
bin, then it would be necessary to determine Mr accurately by considering the i^-correction 
( jWake et al.l l2006l ) . However, our goal here is rather to measure the average Ti2D{f) for 
those galaxies which are brighter than a single magnitude threshold in the redshift range of 
< 2 < 0.4. Hence, a maximum 20% error that may be caused by ignoring iT-correction in 
the estimation of the galaxy absolute magnitude should not be an issue here. 

Fig. [T] plots the observational results of rj2D{f) from the red and the blue galaxy sample 
as solid dots in the top and the bottom panel, respectively. In each panel, the solid line 
represents the analytic model with the best-fit correlation parameters a\ and eni. The errors 
are calculated as standard deviation for the case of no correlation as cr,, = 1 



8np where rZp 

is the number of galaxy pairs belonging to each bin (jPen et al.ll2000l ). As can be seen, clear 
signals of ellipticity correlation higher than Sa^ are detected from both the red and the blue 
galaxy samples at the first bin (r < 3/i~^Mpc). Note also that from the red galaxy sample 
the correlation signal is also detected at the second distance bin (r < 3/;,~^Mpc). 

Table [1] lists the best-fit values of ai and eni. The results show a marked difference 
between the blue and the red galaxy intrinsic alignments: For the blue galaxies, the best-fit 
values of the correlation parameters are found to be ai = 0.20 ± 0.04 and eni ~ 0., which 
indicates that 7720 of the blue galaxies follows a quadratic scaling of ^(r), as predicted by 
the linear tidal torque theory (jPen et al.ll2000l : iLee fc Erdogdul 120071 ). In contrast, for the 
red galaxies, the best-fit values are ai ^ and e^i = (2.6 ± 0.5) x 10~^, which implies that 
ri2D of the bright red galaxies follows a linear scaling of ^(r). This observational result is 
consistent with the picture that the ellipticities of the blue and the red galaxies are induced 
by the initial tidal field and the anisotropic infall of materials, respectively. 
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Before assuring ourselves that the detected signal is real, however, it is necessary to 
examine the effect of the systematic error in the measurement of the position angles of the 
SDSS galaxies because the systematic error could cause correlations of galaxy ellipticities 
that mimic intrinsic alignments. To quantify the systematic error, we perform two simple 
tests. First, we shuffle randomly the redshifts of the selected SDSS galaxies in each sample 
and remeasure r]2D one hundred times. For each distance bin, we calculate the mean averaged 
over these 100 random realizations and the standard deviation between realizations, assuming 
that these 100 shuffling processes are mutually independent. If the systematic error had a 
dominant effect, then a correlation signal would not disappear even when the redshifts are 
shuffled. Fig [1] plots the mean plus and minus one standard deviation as (green) dashed 
lines. As can be seen, when the redshifts of the selected galaxies are shuffled, the ellipticity 
correlations between the galaxies disappear, which suggests that the effect of the systematic 
error on the galaxy intrinsic alignments be negligible. 

Second, we measure the ellipticity cross-correlations between the red and the blue galax- 
ies, which is plotted as solid line and compared with the ellipticity correlation of the red (red 
dotted) and the blue (blue dashed) galaxies in Fig. [21 If the systematic error caused the 
correlation of galaxy ellipticities, then we would find non-negligible cross-correlations of the 
ellipticities between the red and the blue galaxies. As can be seen in Fig. [2], however, 
the cross-correlations are quite weak at the first distance bin, which reassures us that the 
detected ellipticity correlations of the red and the blue galaxies are not spurious intrinsic. 

It is worth discussing the effect of the redshift distortion that we have ignored in our 
analysis. For the blue galaxies, its mean pec uliar velocity disp ersion has been observed to be 



quite small less than 150 km/s on average (jPavis et al.l 119971 ). Therefore it is well justified 



to ignore the redshift distortion effect. On the other hand, for the red galaxies whic h are 



observed to have higher peculiar velocity dispersion around 200 km/s (jPavis et al.ll2003l ). the 
redshift distortion effect may cause non-negligible degree of scatter in the measurement of the 
galaxy separation distance r. Thus, we do not exclude a possibility that the characteristic 
distance scale for the red galaxy intrinsic alignments and the value of ^ni could be smaller in 
real space than observed in redshift space. 

The key implication of our results for the weak lensing survey is also worth discussing 
here. According to our results, the degree of the galaxy intrinsic alignments is quite strong 
for the bright galaxies at redshifts z < 0.4. Since the galaxy intrinsic correlations tend to 



decrease with redshift (ILee fc PenI 120071 ). this implies that at z = 1 relevant for the cosmic 
shear survey the intrinsic correlations may be larger than the detected signal of oi = 0.2 
and thus should not be completely negligible. Our results also suggest that when remov- 
ing close galaxy pairs as intrinsic ellipticity contaminants from the cosmic shear analysis 
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( iHeymans &: Heavenal2003l : iTakada &: Whitell2004j ). a typical distance scale of a galaxy pair 
has to be determined with care, since the scaling of the intrinsic correlations depend on the 
intrinsic property of the galaxies (red or blue). 

We conclude that our first detection of a clear signal of the galaxy intrinsic correla- 
tions from the SDSS data will be useful for the weak lensing measurement as well as for 
understanding the evolution and formation of the red and the blue galaxies. 
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Table 1. The redshift range (z), a total number of galaxies (Ng), the range of the absolute 
magnitude (M^), the best-fit- values of the linear and the nonlinear correlation parameters. 



Color 


z 




Mr 


ai X 10^ 


£nl X 10^ 








(no i^-correlation) 






Red 


[0.0, 0.4] 


283972 


< -19.2 


0.0 ±3.1 


2.6 ±0.5 


Blue 


[0.0, 0.4] 


87188 


< -19.2 


2.0 ±0.4 


0.0 ± 1.6 
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Fig. 1. — Correlations of the two dimensional projected ellipticities from the SDSS red 
and the blue galaxies (top and bottom, respectively) with the absolute magnitude cut of 
Mj. < —19.2 in the redshift range of < 2; < 0.4. In each panel, the solid dots represent 
the observed signals, while the solid line corresponds to the analytic model with the best- 
fit parameters (eq.[2]) found through x^-minimization. The dashed line corresponds to the 
correlation from the data with z-shuffled. For the observed signals, only those galaxies with 
major-to-minor axis ratio oih/a < 0.8 are considered. 
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Fig. 2. — Cross-correlations of the two dimensional projected ellipticities between the SDSS 
red and the blue galaxies (solid). 



